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Abstract 

In this review paper, we consider three kinds of systems of differential 
equations, which are relevant in physics, control theory and other applica- 
tions in engineering and applied mathematics; namely: Hamilton equations, 
singular differential equations, and partial differential equations in field the- 
ories. The geometric structures underlying these systems are presented and 
commented. The main results concerning these structures are stated and 
discussed, as well as their influence on the study of the differential equations 
with which they are related. Furthermore, research to be developed in these 
areas is also commented. 
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1 Introduction 

The aim of this paper is to make evident the relation between some kinds of systems 
of differential equations and certain underlying (hidden) geometric structures, whose 
analysis may allow us to clarify the properties of the system described by those equations. 

We pay attention to three different situations. Two of them concern ordinary differ- 
ential equations (odes), and the third one is partial differential equations (pdes). 

The first one is (autonomous) Hamiltonian systems, with the symplectic geometry 
as background structure. In fact, there is a historical interplay between the equations 
and the geometric structure, and results from one of these aspects depend for their 
interpretation on the other, leading to new insights in both aspects. This subject is 
treated in Sectional 

The second is singular differential equations; that is, those which cannot be written 
in normal form; and in particular, the case where the dependence on the derivatives is 
linear. The underlying structure is a submanifold of a tangent bundle. The study of 
systems of singular differential equations was made separately in theoretical physics and 
in some technical areas such as engineering of electric networks or control theory. In 
both cases, algorithms for solving the problems (which are essentially the same) were 
developed, although independently. The study of the underlying geometric structure 
has proved to be invaluable to find the solutions of these equations and understand their 
properties. All of this is developed in Sectional 

Finally, the third subject presents a different situation, and concerns partial differen- 
tial equations appearing in physics and technical applications. Many partial differential 
equations admit a variational formulation. Multisymplectic geometry underlies this for- 
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mulation. It seems to be a generalization of symplectic geometry, but the difficulty of the 
problems is of increasing magnitude. One of the objectives of this lecture is to present 
some of these problems, enabling us to clarify some aspects related to all these partial 
differential equations. Section HI is devoted to it. 

Manifolds are supposed to be real, paracompact, and C°^. Maps are C°°. Sum over 
crossed repeated indices is understood. 



2 Hamilton equations and symplectic geometry 

2.1 Hamilton equations 

The study of (time- independent) Hamiltonian systems is a classical subject which 
arises from the transformation of the Euler-Lagrange equations, a second-order system 
of ODES, into a first-order system, using the so-called Legendre transformation pQ, [2], 
[100,. 

For a system defined by a Lagrangian function L{t, q\ g*), the Euler-Lagrange equa- 
tions are 

dL d fdL\ 



which is a second-order system for The associated Hamiltonian system 

introduces new variables {pi, . . . ,pm) called momentum coordinates, and the so-called 
Hamiltonian function, defined as H{t, q\pi) = q^Pi — L{t, q\ g*), where (g^, . . . , q"^) are 
supposed to be solved from the relations 

dL . 

With all these data, the Hamilton equations (Lagrange, 1808) are 

, OH OH 
opi dqi 

It is well known that both systems (P) and are equivalent in the following sense: 
if (gi(t),...,g™(t)) is a solution of (P), then {q^{t), . . . ,q'^{t);pi{t), . . . ,p^{t)), with 
dL 

Pi(t) = ■^T-{t,q^(t),q^(t)), is a solution of and conversely. The procedure of trans- 
forming Euler-Lagrange's equations into Hamilton's equations is known as Legendre 
transformation. 

This situation is the usual for classical systems. Nevertheless, there is a large class 
of interesting physical models (mainly related with relativistic systems) for which this 
procedure is not possible. They are called singular systems, and will be discussed in 
Section IHl 

It is important to point out that there are also dynamical systems of Hamiltonian 
type (i.e., whose dynamical equations are of the form (j2I)) which have no Lagrangian 
counterpart; that is, they are not defined by a Lagrangian function, but by giving a 
Hamiltonian function (see, for instance, [Qlj). 
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2.2 Geometric formulation of Hamilton's equations 



At this point we can ask about the geometric structure underlying these systems of 
ODES. For simphcity we can only consider the time-independent case. In the first case, 
we have the tangent bundle TQ of the manifold Q which represents the configuration 
space of the system. The tangent bundle has some natural geometric structures. Among 
them, the vertical endomorphism and the Liouville vector field can be used to construct 
several geometric objects that allow us to express the Euler-Lagrange equations in an 
intrinsic form. Nevertheless, it seems that the geometric structure underlying the Hamil- 
ton equations is more interesting for many applications. In order to obtain it, observe 
that these equations can be written in matrix form 



I \ _ fdH_ dH_ 
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The skew-symmetric matrix lj = ^^-^ ^ ^ can be interpreted as a 2-form u = dg* A dp, 

■ d d 

and the equation means that the contraction of u with the velocity vector q''Tr- +PiT: — 

ag* opi 

OH OH 

equals the 1-form dH = -g—dq^ + -^—dpi. In intrinsic terms, [q^ [t) , pi{t)) is the local 

representation of a path in a manifold M (the phase space of the system), endowed 
with a 2-form u, and C, satisfies the differential equation 

i^u = dH o ^ 

where if: M ^ M is a function. Furthermore, u is required to be symplectic; that is, 
closed and non-degenerate. So we reach the concept of symplectic manifold {M,uj), 
which is a geometric structure similar to that of a Riemannian manifold. In a symplectic 
manifold, given a function / G C°°(M), there exists a unique vector field Xf G X(M) 
associated with /, such that ixf = df. Xf is called the Hamiltonian vector field 
associated with /, and / is called the Hamiltonian function associated with Xf. Hence, 
note that the Hamilton equations for ^ can be written as ^ = X^ o ^, where Xh is the 
Hamiltonian vector field for H. 

Elementary examples of symplectic manifolds are and S"^, both with the area 
element u, and the cotangent bundle T*Q of the configuration manifold Q of the system, 
which carries a canonical symplectic structure. In this last case, the above mentioned 
Legendre transformation is a map J-'L: TQ — > T*Q. 

As we will see in Section 01 this geometric framework can be extended to describe 
singular Hamiltonian systems, just allowing the closed form u to be degenerate (then it 
is called a presjTiiplectic form, and {M,uj) is said to be a presymplectic manifold). 

An skew-symmetric bilinear map on C°°{M), which is called the Poisson bracket, can 
be defined in every symplectic manifold (M, uj) as follows 

{f,g} ■.= ujiXg,Xf) 

which allows us to write the Hamilton equations in the form 

q' = {q\H} , p^ = {p^,H} 
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(Actually, the concept of Poisson bracket is more general than the symplectic structure 

A relevant feature of Hamiltonian vector fields is that they preserve the symplectic 
(resp. presymplectic) form. In fact, X G X(M) is a (locally) Hamiltonian vector field iff 
LxUJ = 0. In the simple example of the symplectic manifold (M^,a;), this result means 
that the area of a region f/ C and the area of the transformed region ip{U) C M^, 
under the flow of a Hamiltonian vector field X, are the same. This flow is an example 
of a map preserving the symplectic form, that is, a symplectomorphism. 

2.3 Main results concerning symplectic geometry 

In 1912, Poincare, working on the symplectic manifold (M^,C(j), proposed a simpli- 
fied 2-dimensional model for the solar system, obtaining the first result on symplectic 
geometry which was used for studying Hamiltonian systems (see j2], P): 

Theorem 1 (Poincare's last geometric theorem). Suppose cp: A ^ A is an area-preserving 
diffeomorphism of the closed annulus A = M./Z x [—1, 1], which preserves the two com- 
ponents of the boundary, and twists them in opposite directions. Then (p has at least two 
fixed points. 

This result leads to the concept of Poincare map for Hamiltonian systems, and to 
determining the existence of periodic orbits for many dynamical systems. There is a 
generalization of this result (the so-called Arnold's conjecture) and several applications 
of it ( 9J, pp. 55-56). In particular, the relation between the number of fixed points of a 
symplectomorphism and the critical points of a Morse function in the manifold [2], |72j . 

This was the standpoint of a long sequence of results on the geometry and topology 
of symplectic manifolds, whose developments have led to a clarification of the structure 
and properties of Hamiltonian systems. Some of them are the following (see, for instance, 
P, [Sn], m, inni for details): 

Theorem 2 (Darboux, 1882). Let {M,uj) a symplectic manifold, and x G M. Then, 
there exists a local chart {U; q^,Pi) at x such that uj\u = dg* A dpi. 

Some consequences of this theorem (which, in fact, predates the Poincare theorem) 

are: 

1. The Hamiltonian systems are locally equivalent in the following sense: given two 
Hamiltonian vector fields, there is a symplectic transformation that maps locally 
one into the other. 

2. There are no local invariants in symplectic geometry (such as curvature in Rie- 
mannian geometry). 

3. The group of diffeomorphisms preserving the symplectic structure is infinite di- 
mensional. 
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An analogous theorem and similar consequences can be stated for presymplectic mani- 
folds. 

Theorems (Marsden-Weinstein's reduction theorem, 1974. jHS])- Let G be a Lie group 
which acts symplectomorphically in a symplectic manifold (M, a;), with associated mo- 
mentum mapping J: M — > g* (where g* denotes the dual of the Lie algebra of G). As J 
is equivariant, ifOEg*isa regular value of J, then J^^{0) is a submanifold of M which 
is invariant under the action of G. Furthermore, if the action of G is free and proper, 
then J~^{0)/G is a symplectic manifold with dimension equal to dim M — 2 dim G. 

The origin of this reduction procedure is very old. The original ideas come from Euler, 
Lagrange, Jacobi, Poisson, Lie and Noether. In modern formulation, Moser, Arnold, 
Guillemin, Sternberg, Marie and many others have contributed to its development and 
applications (see [ZOI and references therein). 

Thus, reduction theory concerns the removal of variables using symmetries and con- 
servation laws and, as a consequence of the theorem, every Hamiltonian system with 
symmetries can be reduced to another Hamiltonian system. In addition, other relevant 
consequences are the Atiyah-GuiUemin-Sternberg and Delzaut theorems 0. 

As above, this theorem can be generalized for presymplectic manifolds, as well as to 
other different singular cases (see, for instance, [22] and the references quoted therein). 

Theorem 4 (Kostant-Kirillov-Souriau, 1970. [ST], |HI])- Let G be a Lie group, andQ its 
Lie algebra. The orbits of the coadjoint action of G in Q* are endowed with a canonical 
symplectic structure u* G f]^{Q*) which is G-invariant. 

This theorem has some important corollaries: 

Corollary 1 If the action of G on (M, cj) is Poissonian (see 15,9f ). then: 

1. The associated momentum mapping J maps the orbits of the action of G in M 
into the orbits j(^) of the coadjoint action of G in Q* . Furthermore, J maps the 
symplectic form uj into the natural symplectic form of . 

2. If the action is also transitive, then J is a local symplectomorphism. 

Roughly speaking, this means that Hamiltonian systems with symmetries are the 
orbits of the coadjoint action. 

Another subject related to this theorem is the so-called geometric quantization: in 
a natural way, a Hilbert space Ti can be associated with every symplectic manifold. 
Therefore, the phase space and the state variable functions are transformed into a Hilbert 
space and self-adjoint operators on it (see, for instance, [SH!, lEHj, [00], |l()2j ). 

Theorems (Arnold-Liouville, 1969. 2J). Let {M,uj, H) be a 2n- dimensional integrable 
system with integrals of motion fi = if, /2, . . . , /„. Let c G M" be a regular value of 
f = (/i, . . . , /„). The corresponding level set /~^(c) is a Lagrangian submanifold of M. 



6 



1. If the flows of Xf-^, . . . ,Xf^ starting at a point p G f^^{c) are complete, then the 
connected component of f^^{c) containing p is a homogeneous space for M."^ . With 
respect to this afflne structure, the component has coordinates (pi, . . . ,(pn, known 
as angle coordinates, in which the flows of the vector fields Xf^, . . . , Xf^ are linear. 

2. There are coordinates ipi, . . . which are known as action coordinates, comple- 
mentary to the angle coordinates, such that the ipi 's are integrals of motion and 
{{pi, . . . , {pn, ipi, . . . , ipn) form a Darboux chart. 

This shows that the dynamics of integrable systems (see jH]) are very simple, and 
they have explicit solutions in these action-angle coordinates (see [Oj, p. 110, and [T^). 

Another interesting result is the classification of normal forms of quadratic Hamilto- 
nians (see p. 486), which is based on WiUiamson's theorem |l()lj : a result about the 
classification of quadratic forms in symplectic vector spaces. 

We also refer to a theorem of Banyaga which states that some classical structures 
and, in particular, the symplectic form, are determined by their automorphism groups, 
the group of symplectomorphisms |¥. Closely related to this result is the Lee Hwa 
Chung theorem which establishes that, in a given symplectic manifold, appart from 
the symplectic and volume forms, the only differential forms invariant by the set of 
locally Hamiltonian vector fields are multiples of exterior powers of the symplectic form. 
The original version of this theorem [44 concerns the uniqueness of invariant integral 
forms (the Poincare-Cartan integral invariant) under canonical transformations (that 
is, those diffeomorphisms in a symplectic manifold, mapping Hamiltonian vector fields 
into themselves), and this result leads to characterize canonical transformations in the 
Hamiltonian formalism of Mechanics, identifying them with the symplectomorphisms. 
Afterwards, these results were generalized to presymplectic Hamiltonian systems |28j . 

There are many other topics concerning symplectic geometry in the realm of Hamilto- 
nian systems. For instance, the study of Lagrangian submanifolds and foliations, which 
has application to several kind of problems 99j; such as, to give a very nice interpretation 
of the Lagrangian and Hamiltonian dynamics and the Legendre transformation [HHl, |94j . 
PH] . to characterize the generating functions of canonical transformations (symplecto- 
morphisms) P], inni, or to get a geometrical framework for the Hamilton- Jacobi theory 
P7j . Apart from this, we can point out the so-called symplectic integrators, which allow 
us to obtain algorithms for the discretization of the Hamilton equations and numerical 
solutions of them, and are based on the conservation of the symplectic form [H], |17j . 
[Snj, inZ!- There is also the study of the Schrodinger equation as an infinite dimensional 
Hamiltonian system [OE] (which is related to geometric quantization). 

Finally, as presymplectic geometry is the arena for singular Hamiltonian systems, 
we must also mention some remarkable results concerning this topic. One of the most 
important is the Weinstein extension theorem [98j, and one of its consequences, the 
coisotropic imbedding theorem [211, ES, which allow to establish the local structure 
of presymplectic Hamiltonian systems as systems defined in a symplectic manifold |14j . 



7 



3 Singular differential equations 



3.1 Systems described by singular differential equations 

In general, an ordinary differential equation is a relation F{t, x, x, x, . . .) =0 involv- 
ing an independent variable t, a dependent variable x(t), and its derivatives up to a 
certain order. When the highest-order derivative can be solved, say x'^'^^ = f(t,x, . . .), 
the equation is said to be in normal form, which is the most appropriate for studying 
and solving the equation, especially thanks to the existence and uniqueness theorem. 
However, we are mainly interested in the case where this highest-order derivative cannot 
be solved, not only on a point, but on an open set, and this is what we will mean with 
the term "singular". 

Although of leeser importance than equations in normal form, singular differential 
equations have been studied for decades, especially in the last 30 years. The main reason 
is that they appear in many applications, under various names: singular, degenerate or 
implicit systems, differential-algebraic equations (DAEs), descriptor systems, . . . Let us 
consider some of these apphcations. 

• Theoretical physics. A Lagrangian L{t, g*, g*) is called singular when its hessian 

d^L _ 

matrix „ . „ . is singular. This means that its Euler-Lagrange equations M are 

singular. Such Lagrangians appear in the description of relativistic phenomena 
and, as a consequence, in all the fundamental theories of physics; they were first 
studied by Dirac and Bergman in 1950. (See, for instance, |12]). 

• Applied mathematics. In control theory, it is usual to deal with linear singular 
systems; as for instance 

Ex = Ax + Bu, y = Cx + Du, 

where E, A, B, C and D are constant matrices, E singular. Such systems appear 
in various problems coming from optimal control, constrained control and electric 
circuits, and so they have been widely studied. One of the first books on these 
systems was [7j; more recently |17j gives a complete account of singular linear 
systems. The nonlinear case has been less studied, and is related to the theory of 
singular perturbations |3B]. 

• Other applications. Many of them can be found in circuit theory, but also in 
chemical and industrial engineering. We also know of applications to more faraway 
fields, such as econometry or biology. 

It is known that the introduction of additional variables allows us to transform any 
differential equation into a first-order autonomous equation, F(x, x) = 0, or x = f(x) 
in normal form;. For the sake of simplicity, from now on we will consider only such 
equations. 
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3.2 Problems arising when solving singular systems 

In general, singular systems cannot be represented by a vector field, so their inte- 
gration — analytical or numerical — leads to new problems, which are not present in the 
regular case. 

Consistency The first point to be noted is that a singular differential equation may 
not have solutions passing through each point of the phase space M of the system. In 
other words, not all values of the variables are admissible initial conditions. 

To solve the equation, first one should identify these admissible initial conditions. 
Usually they are not obtained directly but in an algorithmic way. For instance, the 
differential equation may imply some relations 0"(x) = among the variables x; these 
relations, sometimes called primary constraints, define a subset Mi C M, the primary 
constraint subset. Hence, the first step consists in finding these constraints. 

Nevertheless, the problem does not finish at this point: not only the initial conditions 
Xo must be inside Mi, but also their evolution x(t) must remain there — this is a tangency 
condition. This implies additional constraints, as will be shown later. This procedure is 
repeated until a final constraint subset Mf is obtained. 

Uniqueness Under some favourable conditions of regularity, the preceding procedure 
ends with a submanifold, where the differential equation can be represented by a family of 
tangent vector fields, so that their integral curves describe the solutions. Since different 
vector fields yield different solutions, in general, giving x(0) = Xq, does not determine 
the solution. Sometimes, even the knowledge of all the derivatives at t — does not 
determine the solution. 

This is not always the case. Sometimes the singularity of the initial differential 
equation arises from an "inappropriate" choice of the variables (that is, the initial state 
space is too large to describe the real degrees of freedom of the problem, due to the use 
of redundant coordinates); in such a case, the system has a unique solution on the final 
constraint submanifold. 

Reduction However, in some physical problems the singularity of the system is a con- 
sequence of the existence of a certain kind of internal symmetries (called gauge symmetry 
in theoretical physics), accounting for the fact that there are different solutions repre- 
senting the same physical state. In this case, the differential equation is undetermined: 
for everypoint in the final constraint submanifold which is taken as an initial condition, 
the evolution of the system is not determined because a multiplicity of integral curves 
(of different vector fields solution) pass through it. 

A reduction procedure can be used to remove this ambiguity: different points of Mf 
that can be reached through different solutions beginning at the same initial condition 
must be identified. This quotient is the reduced phase space of the system. Then the 
physical states of the systems are identified with the points of the reduced phase space. 
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Control systems In addition to these problems, other questions are also interesting 
when dealing with singular control systems: controllability, observability, reachability, 
stabilizability, realizability, optimality . . . All these concepts, which are well established 
for regular control systems, must be reconsidered in the singular case. 

Numerical methods Here the matter is the design and convergence of numerical 
methods to solve implicit differential equations, and also their relation with singular 
perturbation problems. There are many articles and some books devoted to these prob- 
lems (see, for instance, jH], [lOj, |iT|). 

3.3 Geometric formulations of singular differential equations 

In geometric terms, a differential equation written in normal form is defined by a 
vector field X on a manifold M. Then the equation for a path 'j: I ^ M reads 7 = X07. 

In the same way, an implicit differential equation can be geometrically described by a 
submanifold D C TM of the tangent bundle of M, and the differential equation is then 
expressed by the inclusion 

7(/) C D. 

Even if the equation is initially set in euclidian space, one is usually led to work on 
submanifolds of the initial space, so the geometric framework is not only nice, but also 
necessary. 

Singular differential equations, from this most general point of view, have often been 
studied in the literature. For instance, in and [HI] this general framework is applied 
to the Hamiltonian dynamics of singular Lagrangian systems. The article studies 
symmetries and constants of motion for these singular equations, whereas their integra- 
bility is studied in [74 . In the same way, but taking M as a euclidean space, [HI] studies 
the existence of solutions of the same problem, and also gives an algorithm for finding 
these solutions under some regularity conditions; this article is indeed a geometric for- 
mulation of the authors' previous works [HD], jEB]- More or less, the same algorithm is 
given in [H2]- 

In the literature, the most widespread singular differential equations are of a special 
type: they combine some restrictions on the base M and some linear relations among 
the velocities in TM. In general, these equations can be written, in coordinates, as 
y4(x)x = 6(x), where A(x) is a matrix, usually singular. Such equations can be called 
linearly singular (they are also called quasilinear) . 

The geometric study of these equations has been developed independently in the 
areas of Theoretical Physics and Applied Mathematics. 

In Theoretical Physics the initial problem was to obtain a Hamiltonian description 
for singular Lagrangians ^H]- The geometrization of this problem was performed later 
|88j . |6Uj . |75j . and a more general framework is that of presymplectic systems [32], where 
the equation of motion has the form 

i^uj = a o 7, 

where is a presymplectic form on a manifold M and a is a 1-form. Hamiltonian 
dynamics corresponds to the case where M C T*Q is a submanifold, u is the puUback 
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of the canonical symplectic form of T*Q to M, and a = dH, where if is a Hamiltonian 
function. An algorithm was presented in the aforementioned paper to study the existence 
of solutions of such a system. Note that, since the transformation X i— > ixuj is not an 
isomorphism, there is no guarantee of either existence or uniqueness of the solutions. 
The same ideas can be applied to search the solutions in the Lagrangian formalism for 
the singular case [30]; but then the presymplectic equation must be supplemented with 
a second-order condition [TH]. Another way to treat this problem is working in the 
manifold M = T*Q XqTQ 

In a more general way, linearly singular differential equations can be described as 
follows a linearly singular system is given by a a vector bundle morphism A: TM — » 
F and a section b of the vector bundle F — > M. Then the differential equation reads 

A o 7 = 6 o 7. 

This problem can also be formulated in terms of vector fields. In this article, Dirac's 
algorithm for presymplectic systems |S2] is generalized to linearly singular systems. This 
framework includes applications to other problems like singular Lagrangian formalism 
and higher order Lagrangians and implicit Hamiltonian systems P|. Symmetries of 
linearly singular systems have been studied recently 36 . 

A less general framework, but still including the presymplectic systems, is presented 
in which corresponds to a linearly singular system with F = T*M; the constraint 
algorithm and symmetries are studied therein. 

Furthermore, some problems in circuit theory led several researchers to consider 
singular differential equations. 

For instance, in [^21 a constrained differential equation is a set of data including 
a linear restriction on the velocities as well as some constraints on the base manifold 
obtained from a certain potential function. The discontinuous solutions of this problem, 
as well as its singularities, are analyzed. 

The same geometric framework for constrained equations is presented in jSH], with 
some simplifications. The problem is to find a curve c: I ^ E such that 

Tp o c = X ° c c{I) C S, 

where p: E ^ B is a. smooth map, x ^ vector field along p, and E is a submanifold of E 
with the same dimension as M. Notice that this problem is also linearly singular. 

We can also mention the paper [16j, which considers a generalized vector field, namely, 
a vector bundle endomorphism A: TM — > TM together with a vector field v in M. The 
equation of motion is then 

y4 O 7 = t> O 7. 

The aim of the paper is to classify the normal forms for generalized vector fields. Later 
papers have studied normal forms [ZB] and stabihty [Hl| of such equations. 

Finally, a brief review on singular system of differential equations is presented in 
where their geometric features, including a geometric method for obtaining a non- 
numerical solution, are analized. 
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3.4 Geometric solution of singular equations 

From a geometric viewpoint, the search for the solutions of an imphcit equation given 
by a submanifold D C TM can be sketched very clearly: since a solution 7 satisfies 
7(i) G -D, necessarily 7(t) belongs to Mi = t{D), where r: TM — ^ M is the natural 
projection. If this is a submanifold (the primary constraint submanifold), we obtain a 
new implicit equation, defined by the subset Di C TMi, where Di = D (1 TMi. This is 
the first step of an algorithm that, assuming appropriate regularity conditions, may lead 
to a consistent differential equation on a certain submanifold Mf. 

The same algorithm applies to a linearly singular system {A: TM ^ F,b), where we 
want to solve the implicit equation A07 = 607. In this case, the equation Ax-X{x) = cr{x) 
for the unknown vector X{x) can be solved only at the points a; G M where this linear 
equation is consistent: b{x) G ImA^;. These points constitute the set Mi. Of course, 
the solutions must be contained in Mi, so it is clear that they are also solutions of the 
linearly singular system (74i:TMi — > Fi,bi) obtained by restricting F, A and b to Mi. 
This is essentially the first step of the constraint algorithm in Dirac's theory. Again, 
under appropriate regularity conditions, this algorithm may lead to a consistent linearly 
differential system on a final constraint submanifold Mf, where bf{x) G ImAf^. Then 
the solutions of the differential equation correspond to the integral curves of a family of 
vector fields Xf + KerAf. Finally, let us remark that this algorithm may be presented 
explicitly in terms of the constraints and vector fields in a suitable way for computation. 

4 Field theory and multisymplectic geometry 

4.1 Geometric formulation of field theory 

As one may see in Tabled many of the main pdes in Physics can be obtained from a 
Lagrangian function. There are other interesting PDEs describing dynamical processes, 
which are also variational with other conditions: constrained problems, higher-order 
Lagrangians, etc. In any case, in all of them, the underlying geometric problems have 
the same level of difficulty as those we are going to analyze next. 

The geometric way of describing these problems, for first-order theories (see, for 
instance, jSHI, 12^1, jSEI, jlH], [HE]), consists in considering a fibered manifold 71: Y ^ 
X (where X is an m-dimensional oriented manifold, with volume form uj G i7"'(X)), 

and the manifold of first-order jets of sections, J^Y ^ Y ^ X. In this situation, 
there are natural geometric structures, such as: the vertical subbundle, the vertical 
endomorphisms, the canonical structure form, the module of total derivations, .... Using 
some of them, we can associate to a Lagrangian function L\ J^Y ^ M a closed form 
VLl G i7™^^(J^y) such that a section ip-.X-^Yis critical for the variational problem 



This equation can be written in other equivalent ways, using m-vector fields or Ehres- 
mann connections [21], [SI]- In a local chart of adapted natural coordinates in J^Y , 




(jV)*iz^^L = 



for every Z G X^J^Y) 
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Name 


Equation 


Lagrangian 


Dirac 


i^^ T / 

7*: Dirac matrices 
I: identity matrix 




Einstein 
(in vacuum) 


Q'. metric ; 7^: Ricci tensor 


i?((?)v/-dct (? 
i?: curvature 


Elasticity 


d'u'^ d fdW\ 

df^ ~ dx^ \du'l J 
H^(x, Vii): stored energy 


1, 




Fluid dynamics 


9u „ 

— + VuU = -gradp 

div u = 
u: fluid velocity; p: pressure 




Klein-Gordon 


V V72 / 2 / 
Qf2 ^ ^ ^ 




(V^)' - m2^2 


Korteweg-de Vries 


du g du d^u ^ 
dt dx dx^ 


ld(j)d(j) 1 
2dx~dt ~ \ 

with 


'd(j)y Id^ 
.dx ) 2 dx'^ 

dd 


Laplace 




i 
2 




Maxwell 


9E 

divE — p , curl H — — — h j 

dt ^ 
dE 

dt 

E TT' plertrir maffnetir fields 

p: electric charge 
j: electric current 


with 

. ^ dA 
E = -grad0 - — 
at 

B = curlA 
0, A: electric, magnetic potentials 


Schrodinger 




2m 

9- 


metric 


Sine-Gordon 


d'i^ . 

= sm 


Idtjj dip 


Wave 






- l|V0f 



Table 1: Some of the main pdes in physics, with their corresponding Lagrangian. 
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(x", y"^, v^), the expression of this condition is the well known system of Euler-Lagrange 
equations for the multiple variational case: 



dL 




d 


/ dL 















= ; A = 1,...,N 



It is interesting to point out that, in the particular case with dimX = 1, this framework 
describes the time-dependent mechanics jSH], IHI], and then the multisymplectic form 
reduces to be a contact or a cosymplectic form. 

As in the case of odes, a Hamiltonian formulation is also possible in this case jlHj . 
[23, 112]; [HH], |ZSI, [HS], after defining the corresponding Legendre map. In all of them, 
the so-called multimomentum phase space is a fiber bundle over Y endowed with a 
multisymplectic structure, which is canonical in some cases, although in others it is 
constructed using additional elements (connections or sections in the bundle). The pdes 
obtained in these formalisms are called the Hamilton-De Donder-Weyl equations. 

To sum up, in all these formulations the underlying geometric structure is a couple 
{M,Q), where M is a differentiable manifold and f2 is a closed /c-form [k > 2), and is 
called a multisymplectic manifold if fl is 1-nondegenerate. 

Natural examples of this structure are and 5" [n > 2) with their volume elements. 
Canonical models of multisymplectic manifolds are a generalization of the cotangent 
bundle of a manifold Q: the multicotangent bundle A^T*Q {I < k < dim Q); that is, 
the bundle of exterior fc-forms on Q, and the subbundles of it (bundles of forms). All 
of them are endowed with a natural multisymplectic {k + l)-form (i.e.; closed and 1- 
nondegenerate) . Other examples are semisimple Lie groups, cosymplectic manifolds and 
Calabi-Yau manifolds (see [45]). 

Finally, it must be remarked that multisymplectic geometry is not in fact the only ge- 
ometrical framework for describing these kinds of variational systems of pdes, and alter- 
native, sometimes equivalent, geometric structures are possible: the so-called polysym- 
plectic jHE], |1H|, k-symplectic jS], [ZH], k-cosymplectic and Lepagean ^3] formalisms. 

4.2 Results and open problems on multisymplectic geometry 

Starting from a multisymplectic manifold {M,Q), one can expect to obtain results 
concerning the geometric structure of M and the pdes defined on M, in an analogous 
way, as in the case of symplectic geometry. Research on these topics has just started 
(see, for instance, jTU], [H], |1S]), and many of those problems are still unsolved. Next 
we review some of these problems and their current status. 

The first fundamental result in symplectic geometry was the Darboux theorem, which 
stated that all symplectic manifolds are locally isomorphic. This result also holds in some 
particular cases of multisymplectic geometry; for instance, when the degree of Q is equal 
to dim M (volume forms). Nevertheless, in the general case, a multisymplectic manifold 
does not admit a system of Darboux coordinates for the multisymplectic form jTlj. In 
fact this is a problem arising from linear algebra: the classification of skew-symmetric 
tensors of degree greater than two is still an open problem. The kind of multisymplectic 
manifolds admitting Darboux coordinates has been identified recently fl], [HZI, and they 
are those being locally multisymplectomorphic to bundles of forms. 
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Nevertheless, further developments have not been achieved. For instance, concerning 
reduction theory, only partial results about reduction by foliations are currently being 
studied The theory of reduction of systems of "multisymplectic" pdes under the 
action of groups of symmetries, obtaining other simpler but also "multisymplectic" ones, 
is under research [12], [70] . 

In the same way, approaches for generalizing symplectic integrators to this geomet- 
ric framework (i.e., the so-called multisymplectic integrators) have only been recently 
developed [U^ . 

Other results similar to those stated in Theorems 0] (Konstant-Kirillov-Souriau) and 
ini (Arnold-Liouville) have not been achieved yet. 

One can expect to see more work on all these subjects in the future. 
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